Quintessence: A Review 

Shinji Tsujikawa 

Department of Physics, Faculty of Science, Tokyo University of Science, 1-3, 
Kagurazaka, Shinjuku-ku, Tokyo 162-8601, Japan 

E-mail: shinji@rs.kagu.tus.ac.jp 

u 

O I 1 Abstract. 

•^ , Quintessence is a canonical scalar field introduced to explain the late-time cosmic 

acceleration. The cosmological dynamics of quintessence is reviewed, paying particular 
attention to the evolution of the dark energy equation of state w. For the field 
potentials having tracking and thawing properties, the evolution of w can be known 
O -1 analytically in terms of a few model parameters. Using the analytic expression of w, 

; I . we constrain quintessence models from the observations of supernovae type la, cosmic 

&Qj microwave background, and baryon acoustic oscillations. The tracking freezing models 

are hardly distinguishable from the A-Cold-Dark-Matter (ACDM) model, whereas in 
thawing models the today's field equation of state is constrained to be wo < —0.7 
(95 %CL). We also derive an analytic formula for the growth rate of matter density 
\q ■ perturbations in dynamical dark energy models, which allows a possibility to put 

Q\ \ further bounds on w from the measurement of redshift-space distortions in the galaxy 

power spectrum. Finally we review particle physics models of quintessence-such as 
those motivated by supersymmetric theories. The field potentials of thawing models 
based on a pseudo-Nambu-Goldstone boson or on extended supcrgravity theories have 
a nice property that a tiny mass of quintessence can be protected against radiative 
corrections. 
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1. Introduction 

The observational discovery of the late-time cosmic acceleration from the Supernovae 
type la (SN la) opened up a new research area in modern cosmology [TJ [2]. About 
70 % of the energy density of the Universe today consists of an unknown component 
called dark energy. This has been also confirmed by other observations- such as Cosmic 
Microwave Background (CMB) [3JH] and Baryon Acoustic Oscillations (BAO) |5j. The 
property of dark energy is characterized by the equation of state w = P/p, where P is 
the pressure and p is the energy density. Dark energy has a negative pressure with w 
less than —1/3. 

One of the simplest candidates of dark energy is the cosmological constant A with 
w = —1. The cosmological constant can arise from a vacuum energy in particle physics, 
but its energy scale is enormously larger than the observed energy scale of dark energy 
[6]. There have been many attempts to construct de Sitter vacua in supersymmetric 
theories. In string theory, for example, huge numbers of de Sitter vacua (~ 10 500 ) can 
be present after the so-called flux compactification of higher- dimensional manifolds [7] . 
We may live in a vacuum with a tiny vacuum energy, but it is generally difficult to 
justify the reason for living such a specific vacuum unless some anthropic principle is 
introduced. So far, it is fair to say that there is no satisfactory scenario where the small 
energy scale of dark energy can be naturally explained by the vacuum energy related to 
particle physics. 

If the cosmological constant problem is solved in a way that it vanishes completely, 
we need to find out an alternative mechanism to explain the origin of dark energy [HJ [9] . 
Broadly speaking, we can classify dark energy models into two classes. The first one is 
based on a specific form of matter- such as quintessence pUl El 121 [131 El E3 HE], k- 
essence [HI [18], and the Chaplygin gas [19]. The second one is based on the modification 
of gravity at large distances (see Refs. [20] for reviews). In both classes the dark energy 
equation of state dynamically changes in time, by which the models can be distinguished 
from the ACDM model. 

Quintessence is described by a canonical scalar field <fi minimally coupled to gravity. 
A slowly varying field along a potential V{4>) can lead to the acceleration of the Universe. 
This mechanism is similar to slow-roll inflation in the early Universe, but the difference 
is that non-relativistic matter (dark matter and baryons) cannot be ignored to discuss 
the dynamics of dark energy correctly. Moreover, the energy scale of the quintessence 
potential needs to be of the order of /?de ~ 10~ 47 GeV 4 today, which is much smaller 
than that of the inflaton potential. 

The dynamics of quintessence in the presence of non-relativistic matter has been 
studied in detail for many different potentials HH [13 HH1 [2U [221 [231 EH [25] . Depending 
on the evolution of w, we can broadly classify quintessence models into two classes [24J: 
(i) thawing models and (ii) freezing models. In the first class, the field is nearly frozen 
by a Hubble friction during the early cosmological epoch and it starts to evolve once the 
field mass drops below the Hubble expansion rate. In the second class, the evolution of 
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the field gradually slows down because the potential tends to be shallow at late times. 
For the inverse power-law potential V(<p) = M i+P cf)~ p (p > 0), there is so-called a tracker 
solution [26] along which w is nearly constant during the matter era and w starts to 
decrease after that. This case belongs to a subclass of freezing models. For thawing and 
tracker models there exist convenient analytic formulas of w J27J EH1 GHJ [30] employed 
to test the models with the data of distance measurements of SN la, CMB, and BAO. 

The redshift-space distortions (RSD) appearing in clustering pattern of galaxies 
[3T| 132] can provide additional constraints on the growth rate of matter perturbations 
8 m . Since the evolution of S m is different depending on the field equation of state 
[331 EU [35] , it is possible to place bounds on w from the data of RSD. In fact there 
exist analytic formulas of 5 m and its growth rate [36], which can be used to constrain 
quintessence models. 

In order to realize the cosmic acceleration today, the mass m^ of quintessence 
(defined by mi = d 2 V (</>) / d(fi 2 ) needs to be extremely small, i.e., \m^\ < H ps lCT 33 eV, 
where Hq is the today's Hubble parameter. In general there is a difficulty to reconcile 
such a ultra light mass with the energy scales appearing in particle physics [37]. 
Moreover, in the absence of some symmetry, the radiative corrections may disrupt the 
flatness of the quintessence potentials required for the cosmic acceleration [38]. However, 
it is not entirely hopeless to construct viable quintessence models in the framework of 
particle physics [39l SQl SH S21 S3! SH SHI SSI S7J . 

In this article, we review several cosmological aspects of quintessence- including its 
cosmological dynamics, analytic solutions of w, observational constraints, and particle 
physics models. The review is organized as follows. In Sec. [2] we present the field 
equations of motion for general quintessence potentials and then proceed to the analysis 
of fixed points for exponential potentials. In Sec. [3] we classify quintessence potentials 
into two classes depending on the evolution of w and then derive analytic solutions of 
w. These solutions are employed to put observational bounds on quintessence at the 
background level. In Sec.H] we derive analytic formulas for the growth rate of matter 
density perturbations and discuss constraints on some of quintessence models from the 
recent data of RSD. In Sec. [5] we review theoretical models of quintessence based on 
supersymmetric theories. Sec. [6] is devoted to conclusions. 

2. Dynamical equations of motion and exponential potentials 

Let us consider quintessence in the presence of non-relativistic matter described by a 
barotropic perfect fluid. The total action is given by 



S = / d 4 x\f^g 



1 -.2 - 1 



-M^R--g^d»<t>dA-V{<P) 



+ S m , (1) 



2 pi 2" 

where g is the determinant of the metric g^ u , M p \ is the reduced Planck mass, R is the 
Ricci scalar, S m is the matter action. We assume that non-relativistic matter does not 
have a direct coupling to the quintessence field <f>. 
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We study the dynamics of quintessence on the flat Friedmann-Lemaitre- Robertson- 
Walker (FLRW) background with the line element ds 2 = —dt 2 + a 2 (t)dx 2 , where a(t) is 
the scale factor with cosmic time t. The pressure and the energy density of quintessence 
are given, respectively, by P^ = (p 2 /2 — V(<f>) and p$ = (p 2 /2 + V((f>), where a dot 
represents a derivative with respect to t. The dark energy equation of state is 

W ^J 2 I 2 ~ V ^. (2) 

P4> 2 /2 + V{(j>) 



The scalar field satisfies the continuity equation p^ + 3H(p ( / ) + P^) = 0, i.e., 

<f> + 3H<j) + Vcp = , (3) 

where H = a/a and V^ = dV/d<p. For a matter fluid with the energy density p m and 
the equation of state w m , the equations of motion following from the action ([T]) are 

3M 2 H 2 = <p 2 /2 + V(<P)+ Pm , (4) 



2M 2 { H 



pi J 

2 



+ (1 +W m )Pm ■ (5) 



In order to deal with the cosmological dynamics of this system, it is convenient to 
introduce the following dimensionless variables [H 



XS * y^VYK (6) 

Vqm p1 h' Vsm p1 h 

The field density parameter Q^ = p^/{3M 2 x H 2 ) can be expressed as 0.$ = x 2 + y 2 . From 
Eq. (J3J) the matter density parameter Q m = p m /(3M 2 l H 2 ) satisfies Q m = 1 — Q^. The 
field equation of state (J2J) reads w = (x 2 — y 2 )/[x 2 + y 2 ). We also define the effective 
equation of state w e s = — 1 — 2H /{3H 2 ), where H/H 2 can be evaluated from Eq. (jSJ) as 

^ = -3x 2 -^(l + Wm )(l-x 2 -y 2 ). (7) 

Taking the derivatives of x and y with respect to iV = In a and using Eqs. (J3J) and (J7J), 
it follows that 

^ = - 3x + ^Ay 2 + ^x [(1 - w m )x 2 + (1 + Wm )(l - y 2 )] , (8) 

J^F = - ^xy + ~y [(1 - w m )x 2 + (1 + w m )(l - y 2 )] , (9) 

where A is defined by A = —M V \V^/V. 

The models with constant A corresponds to the exponential potential [T3| [HI HE] 

V(0) = V e^ /M ^ , (10) 

in which case Eqs. (jSJ) and (jUj) are closed. The fixed points of this system can be derived 
by setting dx/dN = and dy/dN = [T4]: 

• (a) (x, J/) = (0,0), Q^ = 0, w e g = w m , w is undetermined. 

• (b) (x,y) = (±1,0), f^ = 1, w c ff = w = 1. 

• (c) (x, y) = (X/VQ, [1 - X 2 /Q] l/2 ), fy = 1, w cS = w = -1 + A 2 /3. 
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. (d) (x,y) = (^3/2(l+w m )/X 1 [3(l-wl)/2X^) 1 
Q ( p = 3(l + w m )/X 2 , w eS = w = w m . 

If we consider non-relativistic matter (w m = 0), the matter-dominated epoch (w e g ~ 
0, Qj, <C 1) can be realized either by (a) or (d). The point (d) is the so-called scaling 
solution [13J IT4] . along which the ratio fi m /f^ (7^ 0) remains constant. In order to 
realize the matter-dominated epoch by the scaling solution, we require the condition 
A 2 ^> 1. On the other hand, under the condition A 2 < 2, the epoch of cosmic acceleration 
(w e s < —1/3) can be realized by the point (c). This shows that the transition from 
(d) to (c) is not possible, but for A 2 < 2 the system can evolve from (a) to (c). The 
radiation-dominated epoch corresponds to the fixed point (a) with w e s = w m = 1/3. 

In order to study the stabilities of the fixed points (x,y) = (x c ,y c ), we consider 
linear perturbations 5x and 6y about them. Then the perturbations satisfy the following 
differential equations 

£(;)-*(£)■ "-(*£*/*) • (u) 

\ / \ / \ / x=x c ,y=y c 

where fi(x,y) and / 2 (a;,?/) are the r.h.s. of Eqs. (jHJ) and © respectively. If both the 
eigenvalues \i\ and /j, 2 of the matrix M. are negative, the corresponding fixed point is 
stable. If either \i\ or /x 2 is negative, the point corresponds to a saddle. If both \i\ 
and /i 2 are positive, the fixed point is unstable. For complex values of jjl\ and /x 2 with 
negative real parts, the fixed point is called a stable spiral. 

The eigenvalues of the point (c) are //i = (A 2 — 6)/2 and /x 2 = A 2 — 3(1 + w m ) [Ej, 
so that it is stable under the condition A 2 < 3(1 + w m ) for < w m < 1. The condition 
for cosmic acceleration corresponds to A 2 < 2, in which case the point (c) is stable. The 
eigenvalues of the point (a) are fii = — (3/2) (1 — w m ), /i 2 = (3/2) (1 + w m ), so that it is 
a saddle for < w m < 1. This means that, for A 2 < 2, the solution eventually exits the 
point (a) to approach the attractor point (c). 

The dark energy equation of state w for the point (a) is undetermined, but in the 
realistic Universe, x and y are not exactly 0. The early evolution of w depends on the 
initial conditions of x and y. If x 2 ^> y 2 and x 2 ^C y 2 , we have w ~ 1 and w ~ — 1 
respectively. Finally the solution approaches the constant value w = — l + A 2 /3. Since w 
dynamically changes in this way, the quintessence model with the exponential potential 
is observationally distinguishable from the ACDM model. 

For quintessence models in which A is not constant, Eqs. (jHJ) and (JSJ) are not closed. 
In such cases the situation is more involved, but it is possible to derive analytic solutions 
of w by classifying quintessence potentials according to the evolution of w. In the next 
section we shall address this issue. 
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3. Classification of quintessence models and observational constraints 

In order to study the evolution of w for the models with varying A, we derive the 
differential equations for w and Sl^. Using Eqs. ©-Q, we obtain 



w' = (w - 1)[3(1 +w)- Ay / 3(1 + w)fy] , (12) 

% = - 3(tw - w m )n^l - fy) , (13) 

where a prime represents a derivative with respect to iV = In a. Introducing the quantity 
T = VV^/V^, the parameter A obeys 

A' = -^3(1 +w)Q^(T -1)A 2 . (14) 

The evolution of w is different depending on the quintessence potentials and the initial 
conditions. In what follows we discuss three qualitatively different cases: (i) tracking 
freezing models, (iii) scaling freezing models, and (iii) thawing models. In freezing 
models the potential tends to be shallow at late times, which results in the decrease of 
w. In thawing models the mass of the field becomes smaller than H only recently, so 
that the deviation of w from —1 occurs at late times. 

3.1. Tracking freezing models 

For the field density parameter satisfying the relation 

n ij> = 3(l + w)/X 2 , (15) 

w is constant from Eq. ( TT"2l) . If w = w m , then &,</, is constant from Eq. ( fT3l) and hence 
A is constant. This corresponds to the scaling solution (d) discussed in Sec.[2j Since 
the scaling solution is stable for A 2 > 3(1 + w m ) [14], it does not exit to the fixed point 
(c). If A decreases in time, the system can enter the epoch of cosmic acceleration. From 
Eq. f[T4"]) this condition translates into 



r>i. (16) 

The solution ( TT5l) satisfying the condition ( IT6l) is called a tracker [26], along which Vt^ 
increases and hence w < w m . From Eq. ( 1T5|) we have the relation fi'./fl^ = — 2A'/A. 
Using Eqs. (TTB1 and ( TT41) under the condition fl^< 1, the constant equation of state 
along the tracker is [26] 

_ w m - 2(r - 1) 

w = W(o) = ^ZT\ • ( ll > 



For example, let us consider the inverse power-law potential [13j 

V{(j>) = M 4+p <j)- p , (18) 

where M and p (> 0) are constants. Since in this case T = 1 + 1/p > 1, the tracking 
condition ( TT6l) is satisfied. The constant equation of state ( fTTl) is W(o) = (pw m —2)/(p+2) 
and hence w^ = — 2/(p + 2) during the matter era. With the growth of Q^, w starts 
to decrease from w^. Hence the tracker belongs to the class of freezing models. 
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Figure 1. The field equation of state w versus a for the tracker solution. This case 
corresponds to the inverse power-law potential V((f>) — M <f> . The solid curve is 
derived by solving Eqs. (fT2j) and (fl3|) numerically, whereas other curves show the 1-st, 
2-nd, 3-rd order analytic solutions (|21l) . 



The analytic solution f lT7|) is valid in the regime O^ 1. We take into account 
the variation of w by dealing with Q^ as a perturbation to the 0-th order solution (Tl7|) . 
From Eqs. ( 1121) and ( IT4l) . the perturbation 5w around W(q) obeys [30] 

2 d 2 5w 5 — 6u>( ) ddw 9. xr 9 . 2 . . . . . 

a ~da?~ + 2 a ^fa" + 2^ ~ ^(o))^ ~~ 2 W{0) ^ ~ w (P)) n M = ° ' ( 19 > 

where we assumed that T is nearly constant. For Q^(a) we use the 0-th order solution 



QJa) 



Qm\ a 



-3w 



«>) 



ft*o a" 3w(0) + 1 - fi, 



(20) 



where f^ is the today's value of f2^. Substituting Eq. ( 120]) into Eq. f[T9~]) . we obtain the 
following integrated solution 



w a 



w (°) + J2y 



(-ir-^ ( o)(l-<)) 



R/,(a) 



(21) 



;j _ | (n + l)w(o) + 2n(n + 1)W( 0) V 1 ~ ^(«) 

In Fig. [1] we plot the evolution of w derived by the analytic solution (I2TJ) for the 
inverse power-law potential V(4>) = M b cf)~ l . Each curve corresponds to the 1-st, 2- 
nd, 3-rd order solution, whereas the solid curve is derived by solving Eqs. (TT2l) - (TT4l) 
numerically. We find that the analytic solution up to 3-rd order shows good agreement 
with the full numerical result. The analytic expression of w is parametrized by two 
parameters u>( ) and Q^ alone. 

The observational constraints on the tracker model have been carried out in 
Refs. [301 S3 ED]- In addition to the SN la data, the distance measurements of the 



Quintessence: A Review 8 

CMB and BAO peaks provide the information of the background expansion history 
from the recombination epoch to today. From the joint data analysis of Union 2.1 [5T] , 
WMAP7 [52], and BAO (SDSS7 [53J and BOSS [M]), the tracker equation of state 
during the matter era is constrained to be u>( ) < —0.964 (95% CL) under the prior 
W(o) > —1 [50]. For the potential (fTBl this bound translates into p < 0.075. In Ref. [50] 
it was found that the best-fit corresponds to w^ = —1, i.e., the ACDM. If we do not put 
the prior W(o) > —1, the best-fit model parameters are found to be w^ = —1.097 and 
fi^o = 0.717. With the BOSS BAO data [M] the phantom equation of state (w( ) < — 1) 
is particularly favored, but this is not the regime of quintessence. 

3.2. Scaling freezing models 

The scaling solution [131 E] can be regarded as a special case of a tracker along which 
Qtf, = 3(l+u>)/A 2 is constant. During the matter era, w = w m = and hence Q^ = 3/A 2 . 
Since A is constant, T = 1 from Eq. (TT4I) . This case corresponds to the exponential 
potential (TlOl) . but the system does not enter the phase of cosmic acceleration because 
the field equation of state is the same as that of the background fluid. 

This problem can be alleviated by considering the double exponential potential [55] 

V(cj)) = v ie - Xl ^ /M ^ + V 2 e- X2<t>/M * , (22) 

where A, and Vi {i = 1,2) are constants (see Refs. [56j [571 EH] for related potentials). 
For the parameters satisfying the conditions Ai ^> 1 and A2 < 1, the solution first 
enters the scaling regime characterized by f^ = 3(1 + w m )/\\. During the radiation era 
(w m = 1/3) the constraint coming from the big bang nucleosynthesis gives the bound 
Q^, < 0.045 (95% CL) (59J, which translates into the condition Ai > 9.4. The scaling 
matter era (f^ = 3/A 2 , w = 0) is followed by the epoch of cosmic acceleration driven by 
another exponential potential V2e _A2< ^ Mpl . In this case the solution finally approaches 
the fixed point (c) discussed in Sec.[2j 

The onset of the transition from the scaling matter era to the epoch of cosmic 
acceleration depends on the parameters Ai, A2, and V2/V1. The transition redshift is 
not very sensitive to the choice oiV 2 /Vi, so we can set V 2 = V\ without loss of generality. 
In Fig. [2] we show the numerical evolution of w for A2 = with three different values of 
Ai. For larger Ai the transition to w — — 1 occurs earlier. 

The above variation of w can be accommodated by using the parametrization [60] 

. . W„ — Wf /noX 

W{a) =Wf+ l + (a/a t y/r ' ( 23 ) 

where w p and Wf are asymptotic values of w in the past and future respectively, a t is 
the scale factor at the transition, and r describes the transition width (see Refs. [SI] 
for early related works). The scaling solution during the matter-dominated epoch 
corresponds to to p = 0. For A2 = we have Wf = — 1, in which case Eq. (123]) reduces 
to w(a) = — 1 + [1 + (a/a t ) l l T Y l . As we see in Fig. El the parametrization ( 1231) fits the 
numerical solutions of w very well for appropriate choices of a t and r. For the models 
with A2 = the transition width is around r w 0.33, while a t depends on Ai. 
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Figure 2. The field equation of state w versus a for the potential (f22j) with (a) 
Ai = 10, A 2 = 0, (b) Ai = 15, A 2 = 0, and (c) Ai = 30, A 2 = 0. The solid curves 
are the numerically integrated solutions, whereas the dashed curves show the results 
derived from the parametrization (|23|) with w p = and Wf = — 1. Each dashed curve 
corresponds to (a) a t = 0.23, r = 0.33, (b) a t = 0.17, r = 0.33, and (c) a t = 0.11, 
t = 0.32. 



In Ref. [50] the joint data analysis of Union 2.1, WMAP7, and BAO (SDSS7 and 
BOSS) was carried out by fixing r = 0.33. The transition redshift was found to be 
a t < 0.23 (95 % CL). The case (a) shown in Fig. |2] is the marginal one where the model 
is within the 2a observational contour. This shows that w needs to approach —1 in 
the early cosmological epoch. For A2 > the likelihood analysis was also performed 
in Ref. |50j by numerically solving the field equations of motion with suitable initial 
conditions. The model parameters are constrained to be Ai > 11.7, A2 < 0.539, and 
0.256 < Q m o < 0.279 (95%CL). The models with A2 > 0.5 are disfavored because the 
deviation of w from —1 tends to be significant. 

In k-essence models where the Lagrangian P depends on the field <fi and the 
kinetic energy X = —g fJ,u d^(j)d u (j)/2 fT7\ [18] . the condition for the existence of scaling 
solutions restricts the Lagrangian to the form P((f>,X) = Xg(Xe X( ^' Mpl ) [621 163] . where 
g is an arbitrary function in terms of Y = Xe A< ^ //Mpl . The quintessence with the 



exponential potential (P = X 



ce 



-W/M p 



corresponds to the choice g(Y) = 1 — c/Y, 



whereas the choice g{Y) = — 1 + cY gives rise to the dilatonic ghost condensate model 
P = —X + ce x< ^l M ^X 2 [62J. For the multi-field scaling Lagrangian given by 



P( < f> i ,X i ) = Y,Xi9(Xie 



A;</>;/M p p 



(24) 



i=l 



it was shown [M] that a phenomenon called assisted inflation [65] occurs with the 
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effective slope A e ff = (X^=i V-^ 2 ) -1 ^ 2 ; irrespective of the form of g. In the presence of 
multiple fields, the scaling matter era can be followed by the epoch of cosmic acceleration 
even if the individual field is unable to lead to the accelerated expansion. In Refs. [6"6"|l67] 
the cosmological dynamics of assisted dark energy was studied in detail. 



3.3. Thawing models 

In thawing models the field is nearly frozen by the Hubble friction in the early 
cosmological epoch. In this regime one has w ~ — 1, which corresponds to one of 
the fixed points of ffT2|) . The representative model of this class is characterized by the 
potential of the pseudo-Nambu-Goldstone boson (PNGB) [39] : 

U(0) = /x 4 [l + cos(0// a )], (25) 

where /z and f a are constants having a dimension of mass. 

Let us consider the case in which the field initially exists around = 0j and then 
it starts to evolve after the field mass drops below H. We expand the potential V(<fi) 
around — 4>i U P to second order, as V((f>) = J2 n =o V^K&i) (0 ~~ <Pi) n /n\. Using the 
approximation P^ ~ — p^ ~ — V(0j) and redefining the field u = (0 — 0j)a 3//2 , Eq. (J3j) 
reads [281 129] 

1/2 



U — U! U 



-a 3 /V, 



to 



3^(0* 



4 Ml 



-v., 



d 



(26) 



where we assumed 3V (&) / (iM^) > V^,(0j). Provided that \w + 1| -C 1 the evolution 
of the scale factor can be approximated as that of the ACDM model, in which case 



H^a 2 ^}^ + (1 - ^o)a 3 ] from Eq. 



a(t) 



i-a 



1/3 



sin 



h 2 / 3 (t/t 



. Integration of this equation gives 
2M pl 



*A 



Substituting Eq. ( 1271) into Eq. (1261) . we obtain the following solution 

1 — VLm Vrh 



u(t) = Asinh(uit) + -Bcosh(o;£) 



'roO 



t 2 



n 



(60 



uHl - 1 



sinh(t/t/ 



(27) 



(28) 



which is valid for ut\ ^ 1 (i.e., V^,(0j) ^ 0). The integration constants A and B are 
determined by the initial conditions 0(0) = 0.; and 0(0) = 0. Then, the solution is 



0(t) = (f>i + 



V, 



V, 



sinh(o;t) 
tut a sxnh(t/tt 



-1 



Since w + 1 ~ <p 2 (t)/V((pi) under the approximation 



V(<j)i), it follows that 



w + 1 



V, 



ut K V M ((j>i 



2 r 



ut\ cosh(o;t) smh(t/t\) — smh(u)t) cosh(t/tjC) 



-i 2 



(29) 



(30) 



smh 2 (t/t A ) 

The field equation of state can be written as a function of a by using the value wo today 
(t = t , a = 1). Introducing the dimensionless variables 



K = w£a 



3 V& 



F(a) = Jl+[(Sl^)- l -l]a-* 



(31) 
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we obtain the relation cot = K sinh. * ^JcfiVt^/il — VL^q) from Eq. (12"T|) . Then Eq. (1501) 
reads [281 [29] 



w (a) = -l + (l+ Wo ) a 3(x - 1) J r (a) 



where 



J-(a) 



A' 



(K - F(a))(F(a) + 1) K + (K + F(a))(F(a) - 1) 

(k - ^ 1/2 )(^ 1/2 + i)* + (^ + ^o 1/2 )(^ 1/2 - i)* 



(32) 



• (33) 
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Figure 3. The field equation of state w versus a for the potential (f25j) with (a) 
/ /M p i = 0.5, fa/ f a = 0.5 (X - 1.9), (b) f a /M pl = 0.3, &//„ = 0.25 (# = 2.9), 
and (c) fa/Mpi = 0.1, 4 // a = 7.6 x 10~ 4 (_K~ = 8.2). The solid curves correspond 
to numerically integrated solutions, whereas the bald dashed curves show the results 
derived from the analytic solution (J32J with fi^o = 0.73. 



The field equation of state (132]) is expressed in terms of the three parameters w , 
fl^o, and K. The quantity K is related to the field mass squared m| = V^. For the 
potential fl25]) we have K > 1 for < <pi/ f a < n/2 and K < 1 for tt/2 < 4>i/ fa < K, 

[i.e., K 2 < 0), we can derive the similar 

^2^(00/(3^0) -1 PJ- F o r a 



respectively. If 41^^(00/(3^(0*)) > 1 



expression of w by setting X = iK and i\~ = 

phantom field and a thawing k-essence field, analytic solutions of w similar to (|32|) were 

derived in Refs. [2E1 EH] . 

In Fig. [3] we plot numerically integrated solutions of w as well as analytic solutions 
based on ( l32l) for three different values of K. As long as K < 10 and u>o < —0.3, the 
analytic estimation of (|32|) is sufficiently trustable. For larger K, the field mass squared 
|m|| increases and hence the variation of u> around today is more significant. For the 
validity of the Taylor expansion used to derived the analytic solution (I3"2"|) . we require 
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the condition \K — 1| < 0(1). For the models with V^ > 0, the analytic solutions start 
to lose the accuracy for K smaller than 0.5. 

The observational constraints on thawing models have been carried out in Refs. [28| 
|69| [70| [7T| [50] . If the three parameters w , fl^o, and K are varied in the likelihood 
analysis with the prior 0.1 < K < 10, the constraints on i^ are generally weak. After 
the marginalization over K without any prior on wq, Chiba et al. [50] derived the 
bounds -2.18 < w < -0.893 and 0.703 < Q^ < 0.735 (95% CL). If we put the prior 
wo > —1, the field equation of state is constrained to be wq < —0.849 (68 % CL) and 
w < —0.695 (95 % CL). Although there is no statistical evidence that the models with 
Wq > — 1 are favored over the ACDM, the thawing models with — 1 < Wq < —0.7 are 
not yet ruled out observationally. 

4. Constraints from the large-scale structure 

In addition to the background observational constraints discussed in Sec.[3j the 
quintessence models can be distinguished from the ACDM by considering the evolution 
of cosmological perturbations. The peculiar velocity of inward collapse motion of the 
large-scale structure is directly related to the growth rate of the matter density contrast 
5 m . Then the measurement from redshift-space distortions (RSD) of clustering pattern of 
galaxies can constrain the growth history the large-scale structure. The galaxy redshift 
surveys provide bounds on the growth rate f(z) or f(z)a%(z) in terms of the redshift 
z = 1/a— 1, where / = dha.5 m /dha.a and a 8 is the rms amplitude of 5 m at the comoving 
scale 8/i" 1 Mpc (h is the normalized Hubble constant Hq = 100 /ikmsec -1 Mpc -1 ). 
It is then convenient to derive analytic solutions of f(z) and f(z)a 8 (z) to confront 
quintessence models with the observations of RSD. 

Let us consider scalar metric perturbations $ about the flat FLRW background. 
Neglecting the anisotropic stress, the metric in the Newtonian gauge is given by [72] 

ds 2 = -(1 - 2$) dt 2 + a 2 (t)(l + 2<$>)dx 2 . (34) 

We decompose the matter density p m into the background and inhomogeneous parts, as 
Pm = Pm(t) + 5p m (t, x). We also define 

S m = Sp m /p m , 6 m = V 2 v m /(aH) , (35) 

where v m is the rotational-free velocity potential of non-relativistic matter. In the 
Fourier space the matter perturbations satisfy [72j [73] 

S' m = -Z&-9 m , (36) 

H'\ „ / k x :> 



*.+ ^ + 7Tj*" — UJ •■ (37) 

where k is a comoving wave number and a prime represents a derivative with respect to 
N = lna. Taking the iV-derivative of Eq. (13"B"|) and using Eq. (13"T|) . we obtain 



c + ( 2+ §U-(A) 2 $ 



$" + ( 2 + ^ ] $' 



(38) 
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Provided that quintessence does not cluster, we can neglect the contribution of 
quintessence perturbations relative to matter perturbations. In this case the Poisson 
equation is approximately given by 

k \\ 3 



— 1 ^~-Q m S m , (39) 

where S m = 5 m + 3(aH/k) 2 9 rn is the rest-frame gauge-invariant density perturbation. For 
the modes deep inside the Hubble radius {k ^> aH) relevant to the large-scale structure, 
the r.h.s. of Eq. ( 138]) can be neglected relative to the l.h.s. of it. Using the approximate 
relation 5 m ~ 5 m together with Eqs. ([7]) and ( 139|) . we find that Eq. ( |38l) reduces to 

C + \ (1 - 3wfy) S' m - ~n m 6 m ~ . (40) 

In the following we derive analytic solutions for the growth rate / and fa 8 as 
functions of the redshift z. Since Q^ <C 1 at the early cosmological epoch, we expand 
the quintessence equation of state in terms of fl^, as 

oo 

w = w + Y,M%) n - (41) 

n=l 

We introduce the growth index 7, as / = S' m /S m = (ft m ) 7 = (1 -l^) 7 [H E21 E3 • Using 
Eq. (PP with w m = 0, Eq. flUD reads 



3u^(l - fy) ln(l - fy)Jf = ~ - ~w(l - 2 7 )fy + (1 - Q^y - ^(1 - fi^) 1 ^ .(42) 

The solution of Eq. (1421) can be derived by expanding 7 in terms of f^, as 
7 = 7o + Y^n=i ln(Q<t>) n - On using the expansion fj4Tl) as well, we obtain [351 136] 

3(1 -w ) 3 (1-w )(2-3wo) + 2wi(5-6wq) ,. , 2x ,„,x 

7= T^6w7 + 2 (5-6w ) 2 (5-12w ) n * + °W" ( 43 ) 

If w = —1 and wi = 0, then we have 7 ~ 0.545 + 7.29 x lO -3 ^. Since the second 
term is much smaller than the first one, 7 is nearly constant. Even for the models with 
Wq = —1 and W\ = 0.3 (where the value of w today is around —0.8), the variation of 7 
is small: 7 ~ 0.545 + 1.21 x KT 2 ^. 

The relation S' m /5 m = (1 — Q^) 1 can be written in the form 

".^ -<!#. (44) 



Under the approximation that 7 is constant, the term (1 — fi^) 7-1 can be expanded 
around fi^ = 0, as 

(1 - n^y- 1 = i + J2 a ^ n > a - = W - lift - z ) • ( 45 ) 

n=l i=l 

We also expand 1/w in the form 1/w = (l/w )[l + Y^=ifln(Q<t>) n ], where (3 n can be 
expressed in terms of Wj (i = 0, 1, 2, • • •). Then Eq. (1441) is written as 



\llSr 



dflcf, " 3^0^,4, 



i + ^ Cn (^r 



n=l 



n 



^ On-iA , (46) 



i=0 
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with a = (3 = 1. Integration of Eq. ( I4"6"j) leads to the following solution 



) m o exp 



3u>n 



n 






w 



n=l 



(47) 



J mi &S Og 



65 TO , where b is a bias factor. 



where 5 m0 is the today's value of S m . 

The perturbation 5 g of galaxies is related to 5, 
The galaxy power spectrum V g (k) in the redshift space can be expressed as 

V 9 ( k ) = V gg {k) ~ V?V(fc) + ^Vee(k) , (48) 

where \l — k ■ r/(kr) is the cosine of the angle of the momentum vector k to the line 
of sight (vector r), V gg {k) and Veg(k) are the real space power spectra of galaxies and 
9, respectively, and V g g(k) is the cross power spectrum of galaxy-6 1 fluctuations in real 
space. Neglecting the variation of $ in Eq. (|36|) . it follows that 

9 m ~-f5 m . (49) 

The three spectra V gg , V g9 , and V ee depend on (b5 m ) 2 , (bS m )(fS m ), and (fS m ) 2 , 
respectively. Provided that the growth of perturbations is scale- independent, the 
constraints on b5 m and f'5 m at some scale translate into those on ba$ and fas- The 
quantity fa 8 is useful because it does not include the bias factor b. 
Normalizing <5 m0 hi terms of as(z = 0) in Eq. ( 1471) . we obtain [36] 



f(z)<T S {z) = (l-n<y<T S (z = 0)exp 



3w 



In 



n 



00 



a 



+££«« 



("« 



n=l 



The 0-th order solution of fi^ corresponds to w 



Wo, in which case Q 



(0) 



.(50) 



^o(l + 



the first-order solution 



1 + z) 3m, °]. Using the iterative solution w = w + wi^ , we obtain 
JV(1 + ^) 3u,0 [l - ^o + ^o(l + z)* "]"!/^ 



i - n^ + ^ (i + z ) 3w °[i - n q 



n„ 



z \3w ]wi/w Q 



(51) 



We can continue the similar iterative processes, but it is practically sufficient to exploit 
the 1-st order solution (IBTj) for the evaluation of Q^ in Eq. (|50|) . 

There exists a quintessence potential in which w is constant [751 ES] (see also 
Ref. [8]). In this case we have 3 free parameters w , Q^o, a 8 (z = 0) in the expression 
of f(z)a$(z). In tracking quintessence models the coefficients w n (n > 1) are expressed 
in terms of Wo = i0( ), so there are also 3 free parameters wo, Q^o, and o%{z = 0). 
In these cases it was shown that the analytic result (|50|) up to 7-th order terms of c n 
is sufficiently accurate to reproduce full numerical solutions in high precision [36J. In 
thawing quintessence models, when the variation of w is fast at late times, the analytic 
solution ( !50|) is not very accurate unless higher-order terms of c n are taken into account. 

In constant w models the observational data of RSD up to 2012 place the bound 
—1.245 < w < —0.347 (68 % CL), whereas in the tracking models the tracker equation of 
state is constrained to be —1.288 < w^ < —0.214 (68 % CL). Although these constraints 
are still weak, this situation will be improved in future high-precision measurements. 
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5. Particle physics models of quintessence 

There have been many attempts to construct particle physics model of quintessence in 
the framework of supersymmetric theories. Binetruy [76] showed that the inverse power- 
law potential (fl8l) appears in a globally supersymmetric SU(N C ) gauge theory with N c 
colors and the condensation of Nf flavors. In this theory the power p in Eq. ( ITS]) is given 
by p = 2(N C + Nf)/(N C — Nf), which is larger than 2 under the condition N c > Nf > 0. 
Since p is constrained to be smaller than 0.075 [50] . this scenario is not compatible with 
the current observational data. 

In the presence of gravity, any globally supersymmetric theory reduces to a locally 
supersymmetric supergravity theory. In supergravity the four-dimensional effective 
action is given by [77] 



S = / d 4 xy/^g 



-^■R-K^d^d^ifP* -V(<p,<p* 



(52) 



where <p are chiral scalar fields, and K % ^ is an inverse of the derivative of the so-called 
Kahler potential K, i.e., Kij* = d 2 K/d(p i d(p : > . The effective cosmological constant V is 
expressed in terms of K and the superpotential W, as 

V(tp, y?*) = e KIM l [D i W(K ij *)(D j Wy - 3|Wf/M p 2 J , (53) 

where D t W = dW/dip' 1 + {W/M^)dK/dip\ 

The last term in Eq. ( 1531) is negative and hence this can be an obstacle to realize 
a positive vacuum energy required for dark energy. For example, Brax and Martin [78] 
chose a superpotential W = A 3+a (p~ a (motivated by the fermion condensate gauge 
theory mentioned above) and a flat Kahler potential K = (ftp*, but in this case 
the potential V becomes negative for <fi ~ M p \. This problem can be avoided by 
imposing (W) = [7_8j, but such a constraint is generally difficult to be compatible 
with the models of supersymmetry breaking. The Kahler potential of the form 
K = — Mp 1 ln[(<£> + tp*) /M p i], which is present at tree level for both the dilaton and 
moduli fields in string theory, can allow the possibility of canceling the negative term 
(— 3\W p/Mpj). Introducing a new field = (M p i/a/2) ]n((p/M p \) in this case, the kinetic 
term in the action (152]) reduces to the canonical form C^ n = —d fJ- <j)d ll <j)/2. For the choice 
W = A 3+f V" Q , the potential (jHHJ reads [79] 

V{<t>) = MV^W^ 1 , (54) 

where /3 = 2a + 1 and M 4 = M^ 1 A?+ 5 (/3 2 - 3)/2. The positivity of the potential 
requires the condition (3 > ^/3. Then the slope of the exponential potential, A = V^/3, 
satisfies the condition A > a/6. In this case there exists a scaling solution along which 
Q^, = 3(1 + w m )/\ 2 is constant with w = w m , but the potential needs to be modified at 
late times to realize the cosmic acceleration. 

Copeland et al. [79] tried to construct a viable quintessence potential by choosing 

K = M 2 pl [In (<p + pi /M pl ] 2 , W = A 3+a V ~ a , (55) 
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where the field <p is assumed to be real. The kinetic term becomes canonical by 
introducing a new scalar field, <fi = J y/2K ipip * d(p = — (2M p i/3) [1 — \n(2(p / M p \)] ' . 
The field potential is given by 

V = m 4 [2Y 2 + (4a - 1)Y + 2(a - l) 2 ] e ^- Y ) 2 -Mi-Y) ^ Y , (56) 
where m 4 = 2 2a M~ 2 ~ 2a A 6+2a and 

Y = 1 - ln(2^/M pl ) = [-(3/2)(</>/M pl )] 2 / 3 . (57) 

The field exists in the region — oo < < 0, which corresponds to < Y < oo. For \(f>\ -C 
M p i and |0| > M pl , the potential behaves as V oc (-0)" 2/3 and 1/ oc (-0) 2 / 3 e (-*/ M pi) 4/3 
respectively. In the intermediate region there exists a potential minimum with a positive 
energy density. For the initial conditions satisfying \<f>\ <C M p i the quintessence potential 
is approximately given by V((f)) oc (— 0)~ 2//3 in the early cosmological epoch, so that the 
field exhibits a tracking behavior. If the field is initially in the region \<p\ ^> M p \, 
the contribution of the exponential potential is important. In this case a scaling- 
like behavior can be realized during the radiation and matter eras [79]. As the field 
approaches the potential minimum, the Universe enters the epoch of cosmic acceleration. 
A general problem for supersymmetric quintessence models is that supersymmetry 
must be broken if it is to be realized at all in nature. In the gravity and gauge mediated 
scenarios, the supersymmetry breaking is supposed to occur for the energy scale larger 
than (F) 1 / 2 > 10 10 GeV and (F) 1 / 2 > 10 4 GeV (where F 2 is the first term in Eq. ([53]). 
i.e., F 2 = e ' p l D i W{K^*){DjW)*) 1 respectively, to lift the masses of supersymmetric 
scalar particles above 10 2 GeV. In order to give a negligible vacuum energy in Eq. (1531) 
we require that the superpotential takes the form W ~ (F)M P \ ~ m^^M 2 ^ where 
777.3/2 is the gravitino mass (79]. Then the superpotential W = A 3+a (p~ a used above 
gets corrected by the term m^/2M 2 v This gives rise to the correction of the order of 



m 



3 , 2 M 2 j to the quintessence potential, so that the flatness of the potential required for 
the late-time cosmic acceleration can be spoiled. 

Although this problem looks serious, unconventional supersymmetry breaking 
models in string theory may overcome this problem. In Ref. [80] it was suggested that 
we may live in a four-dimensional world with unbroken supersymmetry. In this scenario 
the mass splitting between the superpartners occurs as a result of the excitations of 
the system while maintaining a supersymmetric ground state. Then we do not need to 
worry about the contribution of the supersymmetry breaking terms to the quintessence 
potential. 

There are also some supergravity models in which the above mentioned problem 
can be avoided. In the framework of extended supergravity models [EH E2] the mass 
squared of any light scalar fields can be quantized in unit of squared of the Hubble 
constant H of de Sitter solutions. The de Sitter solutions correspond to the extrema 
of an effective potential V{4>) of a scalar field <fi. Around the extremum at <fi = the 
field potential is given by V((f>) = A + (l/2)m 2 1 2 with A > 0. In extended supergravity 
theories the mass m^ is related to A via the relation m 2 , = nA/(3M 2 j), where n is an 
integer. Since Hq = A/(3M 2 { ) for de Sitter solutions, m 2 = uHq. In the Af = 2 and 
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J\f = 8 extended supergravity theories we have n = 6 and n = — 6 respectively [B2J 1ST] . 
so that the field potentials are 

V(<j>) = 3H 2 M 2 { [1 ± (0/Mpi) 2 ] • (58) 

The energy scale of the supersymmetry breaking is determined by the constant A. If the 
potential ( 158]) is responsible for dark energy, we require that A ps HqM^ ps 10~ 47 GeV 4 . 
The supersymmetry breaking scale is so small that the ultra light mass of the order of 
10~ 33 eV can be protected against quantum corrections. 

The PNGB models based on the potential ( 125]) also allow to protect the light mass 
of quintessence by the U(l) symmetry. An example of a very light PNGB is the so- 
called axion field, which was originally introduced to address the strong CP problem 
[83J. When a global U(l) symmetry is spontaneously broken, the axion appears as an 
angular field with an expectation value (cp) = f a e t( ^^ a of a complex scalar at a scale 
f a . In string theory there are many light axions, possibly populating each decade of 
mass down to the scale H ph 10 -33 eV [SI]. In the limit // — > the potential vanishes, 
so that the symmetry becomes exact. The radiative corrections to V do not give rise 
to an explicit symmetry breaking term because they are proportional to /i 4 . Hence the 
small mass associated with dark energy can be protected against radiative corrections. 

If the PNGB potential ( 125]) is responsible for the cosmic acceleration today, we 
require that Hq ps ji a /Mi\ and hence \i ~ 10~ 3 eV. The field mass squared around = 
can be estimated as m 2 ^ « — {M 2 Jfl)Hl. The slow-roll condition, \M 2 X V^/V\ < 1, 
translates into f a > M p \. Then the field mass is constrained to be \m^\ < H , so that 
the field starts to evolve only recently. As we studied in Sec. 13.31 this belongs to the 
class of thawing quintessence models. 

In supersymmetric theories there have been a number of attempts to explain the 
small energy scale /x ps l(T 3 eV jUJ S21 HSl E3]- Hall et al. [H] tried to relate /i 
with two fundamental scales, the Planck scale M p \ Pd 10 18 GeV and the electroweak 
scale v ~ 10 3 GeV. There is the induced seesaw scale v 2 /M p \ m 10~ 3 eV, which is 
of the same order of \x. If we assume the relation \i m v 2 /M p \ and f a = M p \, 
it follows that |m?| ~ f^/fa ~ t ' 8 /^pi- This gives rise to the mass of the order 
\m^\ w w 4 /M 3 ! « 10~ 33 eV. 

In order to justify the relation \x w v 2 /M p \, Hall et al. [S] proposed supersymmetric 
models with an axion in a hidden sector. In this set up the axion has interactions with 
the quarks q, q c in the form £ int = m q qq c e %< ^^ a at a scale M, where m q is the quark mass 
of the order of the effective supersymmetry breaking scale tub = v 2 /M p \. If at least 
one of the quark flavors has a mass smaller than the order of M, a quark condensate 
forms such that (qq c ) ~ M 3 e 4< ^/ M with an angular field 0. This gives rise to the axion 
potential V = m q M 3 cos(0// a + 0/M), where M is close to tub- Then the scale \i is of 
the order of // ~ m B = v 2 /M p \. 

In summary the thawing models based on the potentials (125]) and ( 158]) are good 
candidates of quintessence from the theoretical point of view. 



Quintessence: A Review 18 

6. Conclusions 

We have reviewed theoretical and observational aspects of quintessence. We classified 
quintessence models in terms of the evolution of the field equation of state w. 

In tracking models the solutions with different initial conditions converge to a 
common trajectory characterized by the analytic solution ( 12~TT) . A typical example of 
this class is the potential V(<p) = M A+P <f)~ p (p > 0), in which case w is nearly constant 
( w (o) = — 2/(p + 2)) during the matter era. The joint data analysis of SN la, CMB, and 
BAO gives the bound u>( ) < —0.964 (95 % CL) and hence the deviation from the ACDM 
is small. The inverse power-law potential appears in a fermion condensate model of a 
globally supersymmetric gauge theory, but the theoretical values of p are larger than 
those constrained by observations. 

The exponential potential V(4>) = Voe~ x ^^ Mpl gives rise to a scaling solution along 
which w = w m and Q^ = 3(1 + w m )/X 2 . Under the condition A 2 > 3(1 + w m ) the scaling 
solution is an attractor during the radiation and matter eras, but it does not exit to the 
epoch of cosmic acceleration. This problem can be alleviated for the double exponential 
potential (|22l) or for the potential ( J56l) appearing in the context of supergravity. The 
likelihood analysis for the potential (1221) with A2 = shows that the transition from 
w = to w = — 1 needs to occur at the early cosmological epoch (a t < 0.23 (95%CL) 
according to the parametrization ( 1231) with w p = and Wf = —1). 

In thawing models there is an analytic solution (1321) of w written in terms of 
the three parameters wq, Q^q, and K. The parameter K is related to the mass of 
quintessence. We require the condition K < 10 to avoid the rapid roll down of the field 
along the potential. Under the prior wq > —1, the today's field equation of state is 
constrained to be wo < —0.695 (95%CL) from the joint data analysis of SN la, CMB, 
and BAO. The potential (1251) of PNGB and the potentials (1581) appearing in extended 
supergravity theories belong to the class of thawing models. In these models, the small 
field mass m^ associated with dark energy can be protected against radiative corrections 
due to underlying symmetries. 

In order to confront quintessence models with the observations of redshift-space 
distortions of clustering pattern of galaxies, we derived analytic formulas for the growth 
rate f(z) as well as f(z)as(z) of matter density perturbations. These are useful to place 
constraints on the quintessence equation of state. We expect that future high-precision 
observations of RSD, combined with other measurements such as SN la, CMB, BAO, 
and weak lensing, will allow us to distinguish quintessence from ACDM. 
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